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Detailed spectrally and spatially resolved radiance has been measured in the Electric Arc 
Shock Tube for conditions relevant to Titan entry, varying atmospheric composition, free- 
stream density (equivalent to altitude) and shock velocity. Permutations in atmospheric 
composition include 1.1%, 2%, 5% and 8.6% CH, by mole with a balance of Nz and 1.5% 
CH. / 0.5% Ar / 98% No by mole, which is consistent with the current understanding of 
Titan’s atmosphere. The effect of gas impurities identified in previous shock tube studies 
were also examined by testing in pure N2 and deliberate addition of air to the CH1/N2 
mixtures. The test campaign measured radiation at velocities from 4.7 km/s to 8 km/s and 
free-stream pressures from 0.1 to 0.47 Torr. These conditions cover a range of potential 
trajectories for flight missions, including a direct ballistic trajectory, a gravity assist fly- 
by, aerocapture or an extremely high-speed entry. Radiances measured in this work are 
substantially larger compared to that reported both in past EAST test campaigns and other 
shock tube facilities. Depending on the metric used for comparison, the discrepancy can be 
as high as an order of magnitude. Due to the difference with previously reported data, a 
substantial effort was undertaken to provide confidence in the new results. With increased 
confidence in the Test 61 data and improved agreement with simulations, potential causes 
for erroneous results in previously report data, such as the effect of oxygen due to Air 
leakage, gas composition and purity, are discussed. The present work provides a new 
benchmark set of data to replace those published in previous studies. The results from 
Test 61 are available at https://data.nasa.gov/ *. 


I. Introduction 


OR certain planetary entry conditions, radiative heating from the bow shock over the vehicle may con- 

tribute sufficient heat load as to influence heat shield design. Codes for predicting radiative heating carry 
significant margin/uncertainty and require validation at conditions relevant to entry. One of the best ways 
to obtain validation data is through testing in shock tubes. The Entry Systems Modeling (ESM) project in 
NASA’s Space Technology Mission Directorate is advancing the state of the art in aerothermal modeling, by 
refining models for radiation and shock layer kinetics. The non-equilibrium radiation from the CN molecule 
is of particular interest to radiative heating during Titan entry. Two previous test campaigns (Test 43 and 
45 conducted from 2003 to 2005') in the Electric Arc Shock Tube (EAST) and two test campaigns in the X2 
shock tube”*® have been shown to provide results inconsistent with simulations. Furthermore, results from 
Test 43 have been shown to be inconsistent with Test 45 and X2 data.” It was reported that the suspect 
results (EAST Test 43) were potentially due to issues with the radiance calibration and carbon contamina- 
tion.?>’ Even though these results have been shown to be potentially questionable, they are still frequently 
used as benchmark data for development of simulation models for Titan entry.” ‘ Since test campaigns 43 
and 45, the EAST facility has undergone substantial upgrades, improving the data quality and quantity 
obtained. With updates to calibration techniques, and improved cleaning of the shock tube reducing carbon 
contamination, it is warranted to conduct tests to update the previously reported data. Due to the limited 
experimental data available prior to the present campaign and the poor level of agreement between previously 
reported data with either Boltzmann, standard QSS approaches or Collisional-Radiative (CR) models,'*:” 
developing a model to accurately simulate non-equilibrium CN radiation has proven to be difficult. 
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II. Description of the EAST Facility 


The EAST facility at NASA Ames Research Center was developed to simulate high-enthalpy, “real 
gas” phenomena encountered by hypersonic vehicles entering planetary atmospheres. In recent years, the 
EAST facility at NASA Ames has been employed for the purpose of obtaining validation data for radiative 
heating. The shock tube produces high velocity flows in gases of known composition and densities relevant to 
various atmospheric entries. These flows are achieved by creating a sudden pressure discontinuity that moves 
hypersonically into the gas in front of it, in the form of a normal shock wave. Because of the short time scales 
involved, the discontinuity does not have time to mix, but rather compresses the gas as it moves forward, 
much like a spacecraft will do to the atmosphere during planetary entry. This shock wave in the driven gas 
of EAST is therefore assumed to be analogous to the bow shock in an entry scenario. The radiating shock 
wave can be imaged as it passes through the shock tube, and separated into different wavelength ranges 
via spectroscopy. The spectroscopic imaging of the shock is important as the radiation varies significantly 
with wavelength. The ability to predict and model both the non-equilibrium and equilibrium radiance is 
dependent upon knowing the spectral features. Having confidence in spectral predictions is important for 
accurate simulations of radiative heating for an atmospheric entry vehicle. 

EAST has the capability of producing super-orbital shock speeds using an electric arc driver with a 
driven tube diameter of 10.16 cm.":!” The region of valid test gas is located between the shock front and the 
contact surface that separates the driver and driven gases. The test duration is defined as the axial distance 
between these two points divided by the local shock velocity. The characteristics of the EAST arc driver 
result in test durations of approximately 4 - 10 ws. Though short, this test duration is often sufficient to 
capture the peak non-equilibrium shock radiation and the decay toward equilibrium conditions. When the 
shocked gas arrives at the location of the test section in the tube, spectrometers attached to Charge Coupled 
Devices (CCDs) are gated, and the spectral and spatial emission of the gas are analyzed. EAST utilizes four 
spectrometers per shot, each associated with four different wavelength ranges. These cameras are referred 
to as: VUV (~ 120 — 420 nm), UV/Vis (~ 190 nm — 500 nm), Vis/NIR (~ 480 nm — 900 nm), and IR 
(~ 700 nm — 1650 nm). For these tests, a small number of measurements were made in the VUV below 190 
nm, but no signal was detected. Therefore these four cameras mostly fixed on four ranges: 190-370 (VUV), 
370-500 (UV/Vis), 500-890 (Vis/NIR), 890-1450 (NIR). The data presented below are therefore considered 
over these ranges, or a subset of these ranges when comparison to previous work is presented. 


II.A. Test Conditions 


The present test campaign conditions were designed to match proposed/previous trajectories for various 
Titan missions and to provide invaluable data for correctly simulating the radiative heat flux and quantifying 
uncertainty margins. Tests were designed to match conditions from the previous (2003 — 2005) EAST 
campaigns, as well as tests conducted in the X2 facility in Australia.? Furthermore, a small number of 
tests also focused on validating higher shocks speeds to open the trajectory space available to future Titan 
mission planners. Figure 1 shows a plot of various proposed and previously flown Titan entry trajectories 
with the Test 61 conditions. These conditions corresponds to a shock speed range of 4.7 to 8 km/s, free 
stream pressures ranging from 0.1 Torr (13.3 Pa) to 0.47 Torr (62.5 Pa) and a test gas composition of either 
95% No / 5% CHa, 98.9% No / 1.1% CHa, 98% No / 2% CHa, 98% No / 1.5% CHy / 0.5% Ar or 100% No 
by mole. The 0.1 Torr condition was the focus of the current work as it has been tested by several other 
facilities and previous EAST campaigns. The 0.47 Torr cases are more likely to correspond to a condition 
close to peak heating for Titan entry. 

At the time of the previous EAST tests (before the Huygens/Cassini mission reached Titan), there was 
a large degree of uncertainty regarding the composition of Titan’s atmosphere, thus testing at both 2% CH, 
and 5% CHy.! Testing at 8.6% CH, was also performed due to a mole versus mass fraction mix up with 
the ordering of the test gas bottle.! Since the entry of the Huygens spacecraft at Titan, it was found that 
the CH, mole fraction was closer to 1.5% with trace amounts of Hy and Ar (< 1% each).'! An approximate 
test gas composition of 1.5 CH4 / 0.5% Ar / 98% No was used to replicate the Titan atmosphere. The 
impact of argon on aeroheating in Titan gas mixtures has not previously been tested experimentally or 
computationally. The radiation specifically from argon is likely to be small, however, argon will increase the 
shock layer temperature, and thus indirectly affect the radiative heating. 
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Figure 1. Comparison of EAST Test 61 testing conditions to various actual and proposed flight missions to Titan. 


III. Summary of Codes 


Experimental data are compared to simulations with DPLR and NEQAIR, used to calculate the flowfield 
and radiation, respectively. This section will provide a brief overview of these codes for reference. 


IIl.A. DPLR 


Deir’? '* uses a finite-volume discretization to solve the reacting Navier-Stokes equations for fluids in 


thermochemical non-equilibrium. While the software was originally designed for steady-state aerothermo- 
dynamic analysis of planetary entry vehicles, DPLR has evolved over the years to include a broad spectrum 
of numerical and physical models that enable it to accurately simulate most compressible flows. Additional 
details on DPLR’s capabilities can be found in the references 12-14. The version of DPLR used in this work 
is v4.04.0. 

A two-dimensional grid of a 3 m sphere with 805 grid points on the stagnation line was used in the 
DPLR solutions. The stagnation line being assumed analogous to the flow behind the shock in the EAST 
experiments. The simulations used a 21-species (CH4, CH3, CH2, HCN, Na, Ni, Cy, He, CH, NH, CN, 
CNT, N, Nt, C, Ct, H, Ht, Ar, Ar*, e~)'° gas model in the computations. In addition to chemical 
non-equilibrium, the flow field was assumed to be in thermal non-equilibrium as well. Consequently, a 
two-temperature (T;,-T-,) model was used. In the two-temperature model employed, the translational and 
rotational modes of the molecules are assumed to be in equilibrium (Ti-=Ttrans=Trot), and are distinct 
from the vibrational and electronic modes of the molecules and the temperature of the free electrons (Te, 
=Ty9=Teiec=Te- ). 


II.B. NEQAIR 


Non-EQuilibrium Alr Radiation (NEQAIR) is a line-by-line radiation code which computes spontaneous 
emission, absorption and stimulated emission due to transitions between various energy states of chemical 
species along a line-of-sight.!° Individual electronic transitions are considered for atoms and molecules, 
with the molecular band systems being resolved for each rovibronic line. Since the report of Whiting et 
al.,!° numerous updates have been incorporated into NEQAIR, such as: using the latest version of the NIST 
atomic database (version 5.0),'’ using the bound-free cross sections from TOPbase,'* incorporating the 
CO» database from CDSD-4000,'” parallelization and improvements to the mechanics of QSS. Of particular 
relevance to this paper was the implementation of a QSS model for CN.?°:?! The current release version of 
NEQaIR is known as v14.0.?? However, due to further improvements to QSS and fixing spurious radiation 
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from atomic nitrogen in the IR, Nz Ist Positive and Nf 1st Negative under certain conditions, the pre-release 
of v15.0 has been used for the calculations presented in this paper. 

An over-prediction by atomic nitrogen in the IR was previously noted to be due to use of the excitation 
rates of Park.” In an attempt to reduce the effect of the over-prediction, a merged database of the rates from 
Huo and Park has been created and discussed by Cruden and Brandis.** The merged database appears to 
have mitigated the issue of providing an over-prediction for atomic transitions. In this work, radiance from 
hydrogen atoms was not simulated. NEQAIR does not have a non-Boltzmann model for hydrogen and the 
radiance predicted for Boltzmann distributions were so large as to provide erroneous integrated radiances. 
As there is no evidence of emission in the EAST spectrum from hydrogen (e.g. from the H-alpha line) except 
at the highest conditions, this approximation is acceptable. 


IV. Nonequilibrium Analysis Methodology 


Insights into the agreement between simulations and experimental results for Earth entry radiation were 
made possible by analyzing integrated equilibrium spectra across a wide range of conditions and conducting 
detailed comparisons of the resulting trends.””’?° For non-equilibrium flows, however, it is not immediately 
clear that one parameter can describe the non-equilibrium intensity observed in several shots over a large 
range of conditions. The temporal intensity changes significantly with shock speed as a result of the un- 
derlying change in the physics. The benefits and drawbacks of three proposed methods for determining a 
parameter that highlights either the relative significance of the non-equilibrium radiative intensity or the 
absolute intensity emitted in the non-equilibrium region were previously analyzed.’ The goal of such a pa- 
rameter is to enable the comparison of several simulations and experimental results. Of these three methods, 
the absolute non-equilibrium metric is used in this work. This metric is computed by integrating the radiance 
within 2 cm of either side of the radiative peak, as shown by the red lines in Fig. 2, and is normalized by 
the shock tube diameter. Computing the metric in this manner has been suggested as a more robust way 
to conduct a comparison as opposed to using the peak intensity, since the comparisons are then not depen- 
dent upon experimental resolution limitations such as gate opening times and spatial smearing due to shock 
movement. Under optically thin conditions, this integral will represent the radiance as observed parallel to 
the direction of the shock. Under optically thick conditions, however, this integral has no additional physical 
meaning and is simply a way of comparing data corresponding to a given optical path-length/shock tube 
diameter. Previous analyses of Titan non-equilibrium radiation were reported using the peak intensity as 
the metric for comparison.”**> As such, some comparisons of the present work to previous data have been 
conducted in this manner. 
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Figure 2. Example of the “Absolute Non-equilibrium Radiance” metric to be used in this work: Integrating the 
intensity from 2 cm before the shock peak until 2 cm after the shock peak. 
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V. Results and Comparisons 


The results of the Titan experiments performed in the recent EAST test campaign will be presented in 
multiple ways, with comparisons against simulations and other experimental data where appropriate. The 
experiments are shown using the non-equilibrium metric across the four spectral ranges that are typically 
used during EAST experiments, as shown in Fig. 3. The data are shown for several different variations in 
methane concentration. Comparisons are also shown with DPLR/NEQAIR simulations, which were performed 
after the test was completed. The level of agreement between simulations and experiments is very good across 
the VUV, UV/Vis and IR spectral ranges, generally within 30%. Even though good agreement is generally 
observed, there is a substantial under-prediction in the Vis/NIR spectral region. It is suspected that the 
reason for the reduced agreement in this range could be due to an issue with the non-Boltzmann modeling of 
CN Red. A less significant under-prediction is also observed in the IR spectral range. The agreement with 
NEQAIR is generally best at the highest (8.6%) and lowest (0%) CH, concentration. The better agreement 
is likely due to a lower uncertainty in the composition of the test gas with either larger CH, percentage or 
pure Nz. Because the test gas is mixed dynamically with mass flow controllers, the uncertainty in the CH, 
fraction is approximately 0.2% regardless of the composition, which is a much larger relative error when the 
CH, fraction is low. The test gas composition is confirmed by conducting Residual Gas Analyzer (RGA) 
readings prior to each experiment. Even though small variations in the test gas composition were observed 
from shot to shot, any variation this caused in the radiation measurements was within the noise of the data. 
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Figure 3. Nonequilibrium Metric Radiance at different compositions of CH,, Nz and Ar at 0.1 Torr. 
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V.A. Comparison with Previous Shock Tube Data 


Results from the recent campaign (Test 61) have been compared with previous Titan experiments performed 
in EAST (Test 43 & 45). Tests 43 & 45 were performed prior to the facility upgrade in 2008. The non- 
equilibrium metric has been used to compare the experimental and simulation results across the different 
campaigns as shown in Figs. 4 to 7. Wavelength regions have been chosen to maximize the number of tests 
than can be used for comparison across tests 43, 45 and 61. Figure 4 shows results for 100% Ng, Fig. 5 shows 
results for 98% Ng / 2% CHa, Fig. 6 shows results for 95% No / 5% CH, and Fig. 7 shows results for 91.4% 
No / 8.6% CHy. Figure 5 highlights that there is a significant discrepancy between Test 61 and Test 45, 
with Test 61 being 2 to 4 times larger in magnitude compared to Test 45. Even though there is a significant 
discrepancy, the trends with velocity between the two data sets are similar. However, as shown in Fig. 7, the 
data from Test 43 shows no discernible trend, therefore comparisons with previous EAST data will focus on 
Test 45. For most velocities and variations of CH, concentration, the simulations agree far better with Test 
61 data. The exception being that for the higher speed shots at 98% Nz / 2% in the 480 to 680 nm range, 
the simulations agree better with Test 45 (see Fig. 5(c)). However, as will be noted later in the paper, the 
spectral comparison is substantially improved when the simulations are compared to Test 61 as opposed to 
Test 45. 
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Figure 4. Non-equilibrium metric comparison at 100% No, 0.1 Torr. 
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Figure 5. Non-equilibrium metric comparison at 2% CH, / 98% No, 0.1 Torr. 


7 of 27 


American Institute of Aeronautics and Astronautics 


. 164 
7) 
“ 
£147 
2 
= 1.2 4 
£ 
c 14 
So 
+ 
TY os 4 
o 
& 
~ 06 4 
2 
s 
® 044 1. 
a = 5 a 
a.) 3 024 . 
a 1 
i= 
6 
T T T T T 1 z 0 T T T T T 
5.5 6 6.5 7 75 8 5 5.5 6 6.5 7 7.5 
Velocity, km/s Velocity, km/s 
(a) 310 — 370 nm (b) 370 — 440 nm 
mTest45 ATest61 -@DPLR/NEQAIR @Test43 mTest45 ATest61 -@DPLR/NEQAIR 
. . és 0.25 
a rf 
A aA Ee 
A A A LL 
0.2 
aA «4 = 
£ 
i= 
o 0.15 
o 
cy 
So 
ive] 
© 01 
2 
a = 
. " 2 
= 0.05 
ov 
a r] 
Mo = 
[) rT] 
7 7 T J T 1 Bw Oo ; r 1 ; : 
5.5 6 6.5 7 7.5 8 5 5.5 6 6.5 7 7.5 
Velocity, km/s Velocity, km/s 
(c) 480 — 680 nm (d) 680 — 860 nm 


Non-eq Metric (310-370 nm), W/cm? sr 


0.35 


0.3 


o 


| Test45 ATest61 -®DPLR/NEQAIR 


5.5 


® Test 43 


6 6.5 7 75 8 
Velocity, km/s 


(a) 310 — 370 nm 


| Test45 ATest61 -®DPLR/NEQAIR 


Non-eq Metric (370-440nm), W/cm? sr 
° b N w 
u Lal ui N ul w uw 


o 


8.5 9 5 


0.35 7 


Non-eq Metric (480-680nm), W/cm? sr 


5.5 


6 6.5 7 75 
Velocity, km/s 


(c) 400 — 430 nm 


® Test 43 


S 

b 

u 
1 


S 
b 
1 


0.1 4 


Non-eq Metric (680-860 nm), W/cm? sr 


= Test45 A Test 61 


6 6.5 
Velocity, km/s 


(e) 680 — 860 nm 


® Test 43 


0.3 4 


= Test45 ATest61 -®DPLR/NEQAIR 


5.5 6 6.5 7 7.5 8 8.5 9 
Velocity, km/s 


(b) 370 — 440 nm 


*Test43 M|Test45 ATest61 -®DPLR/NEQAIR 


5.5 6 6.5 7 7.5 8 8.5 9 
Velocity, km/s 


(a) 480 — 680 nm 


-@ DPLR/NEQAIR 


7.5 8 


Figure 6. Non-equilibrium metric comparison at 5% CH, / 95% No, 0.1 Torr. 
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Figure 7. Non-equilibrium metric comparison at 8.6% CHy / 91.4% No, 0.1 Torr. 
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Figure 8(a) shows the comparison of a radiance profile from Test 61, Test 45 and DPLR/NEQAIR at 
approximately 6.1 km/s, 98% No / 2% CH, and 0.1 Torr in the UV spectral region. The figure shows that 
the DPLR/NEQAIR result compares very well in terms of the peak radiance magnitude with the Test 61 
data. The Test 45 peak radiance is about a factor of 4 lower compared to either Test 61 or DPLR/NEQAIR. 
The figure also shows that there is a significant disagreement between the experimental and simulated decay 
rates, with the simulation showing a slower decay. Figure 8(b) shows a plot of radiance normalized by the 
peak value. Even though the radiance magnitudes for Test 45 and 61 are quite different, the figure shows 
that both tests have very similar decay rates. Figure 8(c) shows a plot of the non-equilibrium metric spectral 
radiance and due to the slower decay rate, the simulated value is larger compared to Test 61. 

The decay rate in a Titan atmosphere is largely controlled by the rate of nitrogen dissociation (Nz + M 
<= N + N+M) coupled with the exchange reaction (No + C == CN + N).*° It was previously reported 
that by increasing the nitrogen dissociation rate improved agreement with the decay rate,”° and this is shown 
to be the case with the new data, see Fig. 9(a). The figure shows that excellent agreement in the decay 
rate is observed when the nitrogen dissociation reaction is increased by a factor of 3 (listed as “N2x3” in 
the figure legend). Increasing the nitrogen dissociation reaction by a factor of 3 is probably at,'° or maybe 
even beyond, the bound of a reasonable uncertainty for this rate. Furthermore, recent quantum chemistry 
calculations suggest that at these temperatures the Park rate provides an over-prediction,”’ indicating that 
some other phenomena, such as the CN non-Boltzmann modeling, might be at play. Figure 9(b) shows the 
impact of applying a Boltzmann state population to the radiation simulation. While increasing the nitrogen 
dissociation rate reduces the peak radiance, the Boltzmann calculation shows that the peak radiance is 
approximately a factor of 3 higher compared to the experiment. This over-prediction by the Boltzmann 
distribution allows for the possibility that the non-Boltzmann model is under-predicting the CN(B) state 
population that is responsible for CN Violet radiation. 
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Figure 8. Simulations compared with Test 61 and Test 45 data at 2% CH, / 98% No, 0.1 Torr and ~6.1 km/s. 
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Figure 9. Variations in simulations to modify peak magnitude and decay rate. 


Figure 10(a) shows a radiance profile comparison of Test 61, Test 45 and DPLR/NEQAIR at approximately 
6.1 km/s, 98% Nz / 2% CH, and 0.1 Torr in the Vis/NIR spectral region. The figure shows that in terms 
of the peak radiance, the DPLR/NEQAIR result is approximately a factor of 1.7 less than Test 61 and 
approximately a factor of 2.5 times larger than Test 45 in terms of the peak radiance magnitude. However, 
unlike the UV spectral region, the figure shows that the experimental and simulated decay rates are very 
well matched. The excellent agreement for the radiance decay rate is highlighted by the normalized radiance 
shown in Fig. 10(b). Figure 10(c) shows the impact of applying a Boltzmann state population to the 
radiation simulation. As with the UV spectral region, the Boltzmann distribution provides a substantial 
over-prediction. This over-prediction by the Boltzmann distribution allows for the possibility that the non- 
Boltzmann model is under-predicting the CN(A) state population that is responsible for CN Red radiation. 
From Fig. 10(d), which shows a plot of the non-equilibrium metric spectral radiance, it is clear that the CN 
Red radiance is under-predicted by the simulation. Even though there is no non-Boltzmann model for C2 
in NEQAIR, good agreement is observed between the Test 61 results and the simulations for Cz Swan (480 - 
520 nm). 

Comparisons with data obtained in the X2 shock tube are compared with EAST results in Fig. 11. Two 
sets of data are shown for X2, and correspond to testing with different shock tube materials and diameters. 
The results labeled “X2 AB” used a steel tube with a diameter of 8.5 cm, while the results labeled “X2 
CJ” used an aluminum tube with a diameter of 15.5 cm. The data are presented in terms of peak radiance, 
as this is the metric applied in the previous studies.2**> In the case of Titan entries, relaxation is slow 
enough in comparison to the shock motion and camera gating time so as to be comparable across different 
facilities, though this is not generally true of all shock tube tests. As can be seen in the figure, results 
from X2 and Test 45 are relatively consistent. When compared to the new EAST data (Test 61), there is a 
significant discrepancy in the peak radiance, up to a factor of 10. This discrepancy was examined in more 
detail, varying many testing parameters and procedures in an attempt to resolve the difference. Factors 
investigated include absolute intensity calibration, methods for filling the test gas into the tube, methods 
for mixing the test gas and the purity of the test gas. Several of these factors are shown in Fig. 11 and are 
discussed in Sect. VI. The level of contamination in each X2 test campaign was reported, and these values 
were used to determine the amount of air added to the Titan test gas in Test 61 for a small number of shots. 
The addition of 1.4% Og significantly reduced the peak radiance, by about 40%, but was not enough to fully 
explain the difference. Further tests were tried with doubling the Oz added to 2.8%, or trying 1.1% CH4 
with 1.4% Og, but these tests still did agree with the X2 or Test 45 results. Figure 11 also shows results from 
Dpir/NEQAIR, and excellent agreement can be seen with the Test 61 data. Results from CR models,~*:*:” 
which were benchmarked to previously reported data, may now provide under-predictions when compared 
to Test 61. 
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Figure 10. Simulations compared with Test 61 and Test 45 data at 2% CHy / 98% No, 0.1 Torr and ~6.1 km/s. 
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Figure 11. Comparison of multiple X2 test campaigns (X2 AB, X2 CJ) with results from EAST tests 45 and 61. 
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V.B. Validation of Calibration with High-Speed Nitrogen 


In order to confirm the calibration procedure with all the optics in place, two high-speed pure nitrogen shots 
were performed. The high-speed condition was chosen to ensure that the shock comes to equilibrium and 
because EAST has previously demonstrated good agreement to equilibrium at these conditions.*” Results 
for the UV/Vis and Vis/NIR spectral regions are shown in Fig. 12. These regions cover the spectral ranges 
of the dominant Titan radiators of CN Violet, CN Red and Cz Swan. The figure shows that very good 
agreement is obtained between the EAST results and NEqQaIR calculations based on CEA*° equilibrium. 
This excellent level of agreement increases confidence in the EAST Test 61 calibration. 
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Figure 12. Radiance measured at 100% Ne, 11.16 km/s, 0.1 Torr compared with equilibrium. 


VI. Examination of Discrepancies 


As was evidently clear in the previous section, there are significant discrepancies between the new ex- 
periments and the previously published data. Therefore, a significant effort was undertaken to investigate 
potential reasons for this discrepancy, with the hope of being able to re-create the old EAST experiments. 
These efforts are broken up into two sections. Section VI.A discusses different modes of facility operation 
and Sect. VI.B discusses contamination and calibration issues. 


VIA. Facility Operation 


Due to the significant discrepancies between the Test 61 and Test 45 data, tests were performed with the 
goal of temporarily stripping back some of the facility upgrades that occurred in 2008. Two of the upgrades 
that were speculated to have an effect on the radiance measurements were the method of gas loading and 
mixing the test gas and the use of an oxygen plasma cleaner. Previous works have shown that the O2 plasma 
cleaner is very successful at reducing carbon contamination in the facility.* 

The current gas loading configuration uses mass flow controls to meter the gas into the facility just prior 
to the shot. This gas loading method allows the gas to be mixed dynamically and ratios to be specified 
with greater flexibility, as opposed to using pre-mixed tanks. Consequently, the delay time between the 
facility loading and shot firing is shortened, which should reduce contamination levels due to outgassing and 
facility leakage. Three variations for filling the test gas were attempted in order to reproduce previous test 
conditions. The first method used a single mass controller with an externally supplied pre-mixed tank of 
98% No / 2% CHy4 by mole (listed as “Pre-mixed” in the figure). The “Pre-mixed” shots would still retain 
a short delay to fire, but may impact the control of the gas ratio as loaded. The second method bypassed 
the gas loading system, and manually loaded the gas directly into the tube (listed as “Manual load” in the 
figures). The “Manual load” shots would have a longer delay time than the first method, using the mass 
flow controller. The third method was the same as the second method, except the manually filled test gas 
was allowed to sit in excess of 10 minutes prior to the shot (listed as “Manual w/Delay” in the figures). 
A final variation combined the manual load process and disabled the O2 plasma cleaner prior to the shot 
(listed as “No O2 plasma” in the figures). This variation would have been the most similar to the practices 
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employed in Test 45, even though it would not replicate the buildup of contaminants over time. Figure 13 
shows the results of these experiments. None of these variations showed a significant enough effect to bring 
the results in line with the Test 45 data. For the most part, data from these tests were within the noise of 
the experiment. Arguably, there might be some reduction in radiance in the Vis/NIR region at low velocities 
when the variants based on manual loading were employed. However, more shots would be needed to verify 
if this phenomenon is real. 
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Figure 13. Non-equilibrium metric comparison at 2% CH, / 98% No, 0.1 Torr with experimental variations. 


VI.B. Contamination and Calibration Issues 


Previous reports have identified potential calibration and contamination issues with pre-upgrade EAST 
results (before 2008). These issues could contribute to the discrepancy. It has been reported from test 
campaigns in other facilities that a small amount of air leakage occurs after isolating the vacuum pumps 
and before firing the test. This amount of air might initially appear small, as it is mostly nitrogen, as is the 
atmosphere of Titan. However the addition of oxygen can be very significant as it will react with the CH, 
to form CO and thereby reduce the carbon available to form CN. To a lesser extent, the oxygen will also 
react with Ny to form NO and thereby reduce the nitrogen available to form CN. Equilibrium simulations of 
these shock conditions show that a small amount of oxygen (1.4%) can reduce the equilibrium CN number 
density by as much as a factor of two. As an example, based on equilibrium calculations, the mole fraction 
of CN is reduced by 40% for a 5.7 km/s and 0.1 Torr condition when 1.4% of Og is added. Therefore, air was 
deliberately added to the EAST test gas to assess this effect. Air was added to provide molar percentages 
of oxygen at 1.4% and 2.8%. The lower of these two levels was chosen to match the reported air leakage 
in the X2 facility,’ and the higher level was selected as an over-test of the air leakage issue. As shown in 
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Sect. V.A, the addition of air resulted in a significant reduction in CN radiance, but not enough to explain 
the discrepancy. 

An additional test examined the possibility of having an incorrectly mixed gas bottle during the campaign. 
While this seems unlikely given that gas bottles are tested and certified by the supplier, there is precedence 
for errors of this nature. In EAST Test 43, a 5% CH, mixture by mass was received from the supplier when 
the tank was expected to be mixed by mole. In Test 61, it was discovered that a pre-mixed tank containing 
a 2% CH, methane mixture had a 0.1% Ar impurity. These type of mistakes would be discovered during 
testing in the current facility operation due to RGA characterization of the test gas performed during every 
shot. However, prior to 2008, no such instrumentation was in use. In the later stage of Test 43 and in Test 
45, the tank was specified as 1.1% CH, mixture by mass, which is equivalent to 2% by mole. A few tests 
were performed using 1.1% CH, (by mole) with a 1.4% Oz impurity added. This test was intended as an 
extreme case of errors in the gas mixture. 

The results of the air addition tests are shown in Fig. 14. The figure shows that a small amount of 
oxygen causes a significant difference, reducing the radiation by almost a factor of 2. The extreme cases of 
2.8% Og or 1.1% CH, with 1.4% O2 produced even lower radiances and are getting close to the reported 
radiation from Test 45. However, even with this unlikely combination of factors, its still does not fully 
explain the discrepancy. The discrepancy can seemingly only be explained by a combination of multiple 
effects, including contaminants in the test gas, radiance calibration issues and inaccuracies in shock speed 
measurements. The accuracy in measuring shock speed in EAST improved after the facility upgrade due to 
more pressure transducers installed in the facility, sampling at a higher rate and the use of a more robust 
technique for calculating the shock speed. The accuracy in shock speed for Tests 43 and 45 was quoted to 
be 1.5%, while the current uncertainty is up to 0.4%.* It was not possible to reproduce the combination of 
effects required to exactly match Test 45 or X2 data during the Test 61 campaign. 
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Figure 14. Non-equilibrium metric comparison at 2% CH, / 98% No, 0.1 Torr with attempts to recreate previous data. 
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VII. Selected Data for Validation 


This section will detail selected experiments from the EAST Test 61 campaign for use in future model 
validation and experimental comparison. The spatially and spectrally resolved radiance is shown for a 
variety of conditions. The 2% CH, / 98% Ng results are shown in this section, with more details and other 
compositions presented in the Appendix. These selected data sets were chosen based on consistency with 
other results from Test 61 combined with good experimental characteristics, such as test time. Figure 15 
shows three experiments at 5.15 km/s, 6.1 km/s and 7.35 km/s across the four spectral ranges measured in 
EAST. Figures 16(a) and 16(b) show the non-equilibrium metric spectral radiance for the UV/Vis and Vis to 
IR spectral regions respectively. Figure 16(c) shows a composite non-equilibrium metric spectral radiance at 
5.15 km/s across all four spectrometers. The relative importance of CN Violet is clearly seen in this figure, 
contributing to over 60% of the total radiance. 
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Figure 15. Radiance measured at 2% CH, / 98% Noe and 0.1 Torr. 


16 of 27 


American Institute of Aeronautics and Astronautics 


Noneq Metric Spectral Radiance, W/cm? um sr 


——T61-20: 7.35 km/s 


—T61-19: 6.1km/s | —T61-16: 5.15 km/s 9, 161-20: 7.35 km/s | —T61-19: 6.1 km/s | —T61-16: 5.15 km/s 


0.1 Torr 
98% N : 2% CH, 


Noneq Metric Spectral Radiance, W/cm? um sr 


300 350 


Figure 16. 


Wavelength, nm 


Noneq Metic Spectral Radiance, W/cm? um sr 


0.1 Torr 
98% N, : 2% CH, 


400 450 500 480 680 880 
(a) (b) 


80 —T61-16: 5.15 km/s 


70 


0.1 Torr 


60 98% N, : 2% CH, 


50 
40 
30 
20 
10 


200 400 600 800 1000 1200 1400 
Velocity, km/s 


(c) 


1080 1280 


Wavelength, nm 


Non-equilibrium metric spectral comparison at 2% CH. / 98% Noe and 0.1 Torr. 


17 of 27 


American Institute of Aeronautics and Astronautics 


VIII. Conclusion 


This paper has detailed the data obtained from the recent Titan test campaign conducted in the EAST 
facility, known as Test 61 and the results are available at https://data.nasa.gov/ >. The results have been 
presented in terms of the peak radiance and non-equilibrium metric as a function of velocity, along with 
profiles of radiance integrated over wavelength and the non-equilibrium metric applied to spectral radiance. 
The data from these experiments have been compared to previously reported data from EAST and two 
campaigns from the X2 facility in Australia. The Test 61 data have been shown to be several factors larger 
in radiance compared to previous experimental data. In some cases, the discrepancy has been shown to 
be up to an order of magnitude depending on the metric used. The experiments have also been compared 
with simulations using DPLR and NEQAIR. Generally, very good agreement has been observed between 
Test 61 and the simulations over most spectral ranges. In particular, the peak radiance of CN Violet and 
the decay rate of CN Red was very well matched. However, the magnitude of CN Red and the decay 
of CN Violet still showed some level of disagreement. Due to the substantial discrepancies observed with 
previously published shock tube data, a substantial effort was undertaken to ascertain what was the root 
cause of the differences. In order to achieve this goal, several variations in the operation of the facility 
were tried. The only option examined that showed improved agreement with the previously published data 
required contamination to be present in the test gas. When air was deliberately added to the test gas 
in Test 61, based on quoted contamination levels for the X2 facility, the CN Violet radiance dropped by 
approximately a factor of 2. However, given all the options tried, it was not possible to fully explain all of the 
discrepancies observed. Finally, experiments showing the best shot characteristics from Test 61 have been 
selected as validation sources for future code validation studies and facility-to-facility comparisons. Due to 
the increased confidence in the Test 61 experiments, it is recommended that the previously reported Titan 
entry data be replaced with the current data set. 
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Appendix: Validation Data 


The appendix presents the radiance profiles and the non-equilibrium metric spectral radiance for the 
experiments selected as validation data. 
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Figure 17. Benchmark data at 100% No2 and 0.1 Torr. 
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Figure 18. Benchmark data at 1.1% CH. / 98.9% Noe and 0.1 Torr. 
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Figure 19. Benchmark data at 2% CHy / 98% Noe and 0.1 Torr. 
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Figure 20. Benchmark data at 2% CHy / 98% No and 0.1 Torr. 
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Figure 21. Benchmark data at 2% CHy / 98% Noe and 0.1 Torr. 
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Figure 22. Benchmark data at 1.5% CH / 98.5% Nz / 0.5% Ar and 0.1 Torr. 
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Figure 23. Benchmark data at 5% CHy / 95% Noe and 0.1 Torr. 
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Figure 24. Benchmark data at 8.6% CH. / 91.4% Noe and 0.1 Torr. 
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